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Abstract 

Using a sample of (225.2±2.8) x 10 6 J/ip events collected with the Beijing Spectrometer (BESIII) 
at the Beijing Electron-Positron Collider, CP and P violating decays of rj, rf and rj c into ir + ir~ 
and 7r 7r° are searched for in J/ip radiative decays. No significant ij, rf or r\ c signal is observed, 
and 90% confidence level upper limits of B(r) -> 7r + 7r~) < 3.9 x 10~ 4 , B{rf ->■ tt + tt^) < 5.5 x 10~ 5 , 
B(i] c -> vr+vr^) < 1.3 x 10~ 4 , B(r} -»■ vr°7r°) < 6.9 x 10" 4 , B(rf -> 7r°7r°) < 4.5 x 10~ 4 and 
B(rj c ->■ 7T 7r°) < 4.2 x 10" 5 are obtained. 

PACS numbers: 11.30.Er, 13.25.Gv, 14.40.Be 
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I. INTRODUCTION 



Finding the source of CP violation is one of the most important goals of particle physics. 
Violation of CP symmetry has important consequences; it is one of the key ingredients 
for the matter- antimatter asymmetry in our universe. CP violation can be experimentally 
searched for in processes such as meson decays. The decays T]/t]'/t]c vr + 7r~ and 7r°7r°, 
which violate both P and CP invariance, provide an excellent laboratory for testing the 
validity of symmetries of the physical world. In the Standard Model (SM), such decays can 
proceed only via the weak interaction with a branching fraction of order 10 -27 according to 
Ref. Higher branching fractions are possible either by introducing a CP violating term 
in the QCD lagrangian (a branching fraction up to 10 -17 can be obtained in this scheme) 
or allowing CP violation in the extended Higgs sector (in this case 10" 15 can be reached), 
as described in Ref. jl|. The detection of these decays at any level accessible today would 
signal P and CP violations from new sources, beyond any considered extension of the SM. 

The best previously published results for the i], rj' and rh. decays to tt + 7t~ and 7r°7r° are: 
B(r] -»■ 7T+7T-) < 1.3 xl(T 5 Q, B(rf tt+tt") < 2.9 xlO" 3 B{rj c -»■ ti + tj-) < 6xl0" 4 tl, 
B(rj ->■ 7r°7r°) < 3.5 x 10" 4 [J, B[rf ->■ 7r°7r°) < 9 x 10" 4 [g], and B(rj c -> vr 7r°) < 4 x 10" 4 4| 
at the 90% confidence level (C.L.), respectively. 

In this article, results are presented on direct searches for the decays of f]/f]' /rjr — > 7r + 7r~ 
and 7r°7r° with the BESIII experiment based on (225.2 ± 2.8) x 10 6 J /if) events [7|. 



II. BESIII AND BEPCII 

n 

BESIII/BEPCII [8] is a major upgrade of the BESII experiment at the BEPC accel- 
erator j^J for studies of hadron spectroscopy and r-charm physics [loj]. The design peak 
luminosity of the double-ring e + e~ collider, BEPCII, is 10 33 cm _2 s _1 at a beam current of 
0.93 A. The BESIII detector with a geometrical acceptance of 93% of 4ir, consists of the 
following main components: 1) a small-celled, helium-based main draft chamber (MDC) 
with 43 layers. The average single wire resolution is 135 /im, and the momentum resolu- 
tion for 1 GeV/c charged particles in a 1 T magnetic field is 0.5%; 2) an electromagnetic 
calorimeter (EMC) made of 6240 Csl (Tl) crystals arranged in a cylindrical shape (barrel) 
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plus two endcaps. For 1.0 GeV photons, the energy resolution is 2.5% in the barrel and 
5% in the endcaps, and the position resolution is 6 mm in the barrel and 9 mm in the 
endcaps; 3) a Time-Of-Flight system (TOF) for particle identification composed of a barrel 
part made of two layers with 88 pieces of 5 cm thick, 2.4 m long plastic scintillators in each 
layer, and two endcaps with 96 fan-shaped, 5 cm thick, plastic scintillators in each endcap. 
The time resolution is 80 ps in the barrel, and 110 ps in the endcaps, corresponding to a 
K/tt separation by more than 2a for momenta below about 1 GeV/c; 4) a muon chamber 
system (MUC) made of 1000 m 2 of Resistive Plate Chambers (RPC) arranged in 9 layers 
in the barrel and 8 layers in the endcaps and incorporated in the return iron yoke of the 
superconducting magnet. The position resolution is about 2 cm. 

The estimation of physics backgrounds is performed with Monte Carlo (MC) simulations. 



The GEANT4-based simulation software BOOST [llj includes the geometric and material 
description of the BESIII detectors and the detector response and digitization models, as 
well as the tracking of the detector running conditions and performance. The production 
of the J/ip resonance is simulated by the MC event generator KKMC 12], while the decays 
are generate d by EvtGen 13) for known decay modes with branching fractions being set 
to the PDG [y| world average values, and by Lundcharm Ijj] for the remaining unknown 



decays. The analysis is performed in the framework of the BESIII Offline Software System 
(BOSS) [3] which takes care of the detector calibration, event reconstruction and data 
storage. 



III. SEARCH FOR rj, rj' AND rj c DECAYS INTO tt+tt" 



To search for rj, rj' and r] c decays into 7r + 7r~ in J/ip radiative decays, candidate events 
with the topology 77r + 7r~ are selected using the following criteria. Charged tracks are 
reconstructed from MDC hits. To optimize the momentum measurement, we select tracks 
in the polar angle range | cos^l < 0.93 and require that they pass within ±10 cm of the 
interaction point in the beam direction and within ±1 cm of the beamline in the plane 
perpendicular to the beam. Electromagnetic showers are reconstructed by clustering EMC 
crystal energies. Efficiency and energy resolution are improved by including energy deposits 
in nearby TOF counters. Showers identified as photon candidates must satisfy fiducial and 
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shower-quality requirements. The minimum energy is 25 MeV for barrel showers (| cos#| < 
0.8) and 50 MeV for endcap showers (0.86 < |cos#| < 0.92). To exclude showers from 
charged particles, a photon must be separated by at least 20° from any charged track. EMC 
cluster timing requirements suppress electronic noise and energy deposits unrelated to the 
event. 

The TOF (both Endcap and Barrel) and dE/dx measurements for each charged track 
are used to calculate Xpid(^) values and the corresponding confidence levels Probpmii) for 
the hypotheses that a track is a pion, kaon, or proton, where % [i = n/K/p) is the particle 
type. For pion candidates, we require ProbpjD^) > ProbpjD(K) and Probpioin) > 0.001. 

Candidate events must have two charged tracks with zero net charge, and the number of 
photons should be greater than or equal to one. The photon candidate with the maximum 
energy in the e + e~ center-of-mass (CM.) frame is taken to be the J/ip radiative decay photon 
('jrad)- At least one charged track must be identified as a pion. We do a four-constraint 
(4C) kinematic fit imposing energy and momentum conservation under the J/i[) — > 77r + 7r~ 
and J/ip >yK + K- hypotheses and require X^ +7r - < 30 and X^+ n - < X^k+k-- 

Due to the large branching fraction of J/ip — > pir, the decay channel J/ip — >■ pir — > it + it~it 
constitutes the main source of background. To suppress it and other possible backgrounds 
with a 7T°, for candidate events with two or more photons, all pairings of the radiative photon 
and the remaining photons in the event are used to form possible 7r° candidates, and the 
pairing with its mass closest to the n° nominal mass is selected. A clear n° signal is observed. 
To remove it, rn llrad < 0.12 GeV/c 2 or m llrad > 0.15 GeV/c 2 is required. The efficiencies of 
this requirement are 99% for 7](t]') — > it + it~ and 96% for r\ c — > it + it~. 

To suppress the J/ijj — > e + e~ background and other possible backgrounds with e + e~, 
the requirements of E^ MC < 1.2 GeV and \XdE/dx{^)\ < 3 are imposed for both n + and 
7r~ candidates. Here E^ MC means the deposited energy of pion candidates in the EMC. 
To suppress J /if) -> // + //~ background, we require E^ c > 0.4 GeV or E^ MC > 0.4 GeV 
only in i] c — > 7r + 7r~ since its contamination in the low 7r + 7r~ mass region is very small. This 
requirement removes 99.96% of the J/vp — > fi + fi~ background events, while the efficiency is 
62.4% for i] c — > 7r + 7r~ according to MC simulations. 

To avoid possible bias, we first did the background analyses with the signal regions (m^ = 
0.53-0.56 GeV/c 2 for r] -> tt+tt", 0.95-0.97 GeV/c 2 for rj 1 ->■ tt+tt- and 2.95-3.02 GeV/c 2 for 
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r] c — > 7T + 7r ) blinded to check if all the simulated backgrounds described the data outside of 
the signal regions. The same was done in the analysis of 77, rf and rj c decay into 7r°7r° (m n7T = 
0.52-0.57 GeV/c 2 for 77 ->• 7r°7r°, 0.93-0.98 GeV/c 2 for 77' ->• ttV and 2.95-3.02 GeV/c 2 for 
i] c — y 7T°7r ). Backgrounds from a number of potential background channels listed in the 
PDG [lij are studied with MC simulations. An inclusive J/tp MC event sample is also used 
to investigate other possible surviving background events. 

The main non-peaking backgrounds are from J/ip — > pir, J/ip — )■ J/ip — >■ e + e~, 

J/ip a 2 (1320)7r 77T+7T-, J/ip &i(1235)tt ^tv+tv-, J/^j tt+tt-, J/ip 
7o-// 2 (1270)//o(1500)// (1710) ->■ 7?t + 7t-, and ISR events e+e" -> *y ISI £K + ir . The pos- 
sible peaking backgrounds are J/ip —> '-ft] with 77 — >• 7?t + 7t~ for 77 — > tt + tt~, — >• 777' with 
77' — >■ 7p° — > 77r + 7r _ for 77' — > 7t + tt~, and J/if) 777,, with 77 c — )■ 7?t + 7t~ for i] c — > 7t + tt~. 
Branc hing fractions of J ftp — > 777 — >• 77tt + tt~ and J/ip 77/ — > 77?t + 7t~ have been mea- 
sured 14j. After event selection, the contribution from the decay J/ip — > 777 — > 77?T + 7r~ 
to the background is less than 11 events within a mass region of ±2(7 around the 77 mass 
peak. Since the tt + tt~ mass distribution is smooth, this background may be neglected. The 
contamination from the J /if) — > 777' — > 77?t + 7t~ background channel will be fixed in the fit 
below. As for the possible peaking background r] c — > 7tt + tt - for 7] c — > tt + 7t~, this process 
is OZI suppressed, unlike 77/77' —t- 7?t + 7t~, and can only happen by cc annihilation, and the 
photon must be from the final state quark. The branching fraction for such a process is 



expected to be very small. It was calculated to be B(rj c — > 7?T + 7r~) = 4.5 x 10 -6 in Ref. 17] 
in the framework of nonrelativistic QCD. After QED contributions are taken into account, 
the £>(77 c 7tt + 7t~) becomes 1.5 x 10 -6 . With these calculated branching fractions, the 
contribution from this peaking background after event selection can be neglected. 

Figure [T] shows the 7r + 7T _ invariant mass distributions of the final candidate events in the 
77, 77' and 77c mass regions after removal of the blinded boxes, where dots with error bars are 
data and the dashed histograms are all the simulated normalized backgrounds. No evident 
signal is observed, and the simulated backgrounds describe the data well. For 77 — > ty + ti~ , a 
fit with an 77 signal shape obtained from MC simulation and a 2nd order Chebychev function 
as the background shape gives 17 ± 23 signal events with a statistical significance of 0.8<r. 
For 77' — > 7r + 7r~, by fitting the tt + 7t~ invariant mass spectrum with the MC determined 
shape for the 77' signal, the normalized mass distribution from J/tp — > 777' — >• 77tt + tt - as 
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the only peaking background, and a 2nd order Chebychev function for other backgrounds, 
0.1 ± 15 events are obtained with a statistical significance of O.lcr. For rj c — > 7r + 7r~, a fit 
with an acceptance-corrected r) c signal shape obtained from MC simulation and a 3rd order 
Chebychev function as the background shape gives 52 ± 35 signal events with a statistical 
significance of 1.5er. The fitted results are shown in Fig. [T] with solid lines, where the arrows 
show the signal mass regions which contain 95% of the signal according to MC simulations. 
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FIG. 1: (a)-(c): The ir + ir~ invariant mass distributions of the final candidate events in the r/, 
7/ and r\ c mass regions, respectively. The dots with error bars are data, the solid lines are the 
fit described in the text, and the dashed histograms are the sum of all the simulated normalized 
backgrounds. The arrows show mass regions which contain around 95% of the signal according to 
MC simulations. 



We determine a Bayesian 90% confidence level upper limit on N S i g by finding the value 



Mf such that 



Io NsiB CdN t 



sty 



0.90, 



fa £dN sig 

where N sig is the number of signal events, and £ is the value of the likelihood as a function 
of N S ig. The upper limits on the numbers of rj, 77' and r\ c are determined to be 48, 32, and 
92, respectively. 



IV. SEARCH FOR rj, rf AND r) c DECAYS INTO 7r°7r° 

To search for 77, 7/ and r/ c decays into 7r°7r° in J/ip radiative decays, candidate events 
with the topology 77r°7r° are selected using the following selection criteria. An event must 
have 5 or 6 photons and no charged tracks. Here in selecting photons, the EMC cluster 
timing requirement is not used. For the rj/i]' — > 7r°7r° modes, the photon with the maximum 
energy is identified as the radiative photon, and all the other remaining photons in the 
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event are used to form two tt° candidates. For r] c — > 7r°7r°, for photons with E < 0.3 GeV 
(a potential radiative photon), all possible two photon pairings of the remaining photons in 
the event are used to form two tt° candidates. If there is more than one radiative photon 
candidate (E < 0.3 GeV), the one that gives the smallest |m 57 — mj/f\ is used. The 
candidate event is chosen as the photon pairing combination giving the minimum \ = 
a/ (m 77 i — m^o) 2 + (m 77 2 — m T o) 2 . Endcap-endcap combinations are of lower quality and 
only make up about 0.7% of all the tt° candidates, so they are removed. The 77 invariant 
masses are required to satisfy |m 77 — 7?V)| < 0.01625 GeV/c 2 (~ 2.5a) for 77(77') — > tt 7t° and 
|m 77 — 77V) I < 0.0175 GeV/c 2 (~ 2.5a) for r] c — > 7r°7r°. To improve the 7r°7r° mass resolution, 
especially for the 77(77') — > 7r°7r° mode, a 4C kinematic fit under the J/ip — > 77r°7r° hypothesis 
is done. To suppress backgrounds with a u(— > 771°), events with the invariant mass of 771° 
satisfying 0.72 < m 77r o < 0.82 GeV/c 2 are rejected. 

The main backgrounds with tt° signals are from J/ip — > a2(1320)?T — > 77T°7r°, 
J /if; -> 6i(1235)tt ->• 77T°7T°, J /if; W7T°(7T ) with u ->■ 77T°, and J/> ->■ 

7tT// 2 (1270)//o(1500)//o(1710)// 4 (2050) ^ 7tt 7t°. Non-7T° background contributions are 
represented by the normalized number of events in the sidebands of the two tt° mass distri- 
butions, where the two 7r° mass sidebands are defined as 0.084 < m 77 < 0.098 GeV/c 2 or 
0.168 < m 77 < 0.182 GeV/c 2 . 

Figure [2] shows the 7r°7r° invariant mass distributions of the final candidate events in the 
77, 7/ and 77c mass regions with the blinded boxes removed, where dots with error bars are 
data and the dashed histograms are all the simulated normalized backgrounds. No evident 
signal is observed and the simulated backgrounds describe the data well. Since there are 
no peaking backgrounds in r]/r]'/r] c — > 7t°7T°, the invariant mass distributions are fit with 77, 
77' and acceptance-corrected rj c signal shapes obtained from MC simulations and 2nd order 
Chebychev functions as background shapes, and 11 ± 18, 75 ± 30 and 0.1 ± 14 rj, rj' and Tfc 
signal events, respectively, are obtained with the corresponding statistical significances of 
0.6cr, 2.6cr and O.lcr. The fitted results are shown in Fig. [2] as solid lines. 

In the 77' — > 7r°7r° invariant mass distribution, a few events accumulate in the signal region. 
They may from / (980) decays since the process J/ip —> 7/o(980) — > 77t tt is not forbidden, 
and the /o(980) mass is not very far from the small peak. But since the signal significance 
is only 2.6cx, and in the PWA of J/ip — > 77t tt / (980), no signal is evident or needed 18J, 
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the peak may be due to a statistical fluctuation. To give a conservative measurement on the 
upper limit and since this was a blind analysis with the background estimated before hand, 
the process J '/if) — >■ 7/0 (980) — > 77r°7r° is not considered in the fit. 
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FIG. 2: (a)-(c): The 7r°7r° invariant mass distributions of the final candidate events in the r], rj 
and rj c mass regions, respectively. The dots with error bars are data, the solid lines are the fit 
described in the text, and the dashed histograms are the sum of all the simulated normalized 
backgrounds. The arrows show mass regions which contain around 95% of the signal according to 
MC simulations. 



V. SYSTEMATIC UNCERTAINTIES 

The sources of the systematic errors for the branching fraction measurements are summa- 
rized in Table [B The uncertainty is negligible for pion identification since the identification 
of onl y o ne of the pions is required. The uncertainty in the tracking efficiency is 2% per 



track [l9j, and is additive. The uncertainty associated with the kinematic fit is determined 
to be 2% for T]/r]'/r] c — > 7r + 7r~ using the control sample J /if) — > 7r + 7r~7r°, while the uncer- 
tainty can be neglected for r)/r)'/r] c — > 7r°7r° since we require only a very loose \ 2 requirement 
(the default value is 200). The uncertainty due to photon detection is 1% per photon. This 
is determined from studies of photon detection efficiencies in well understood decays such as 
J/ifj p°7T and a study of photon conversion via e + e _ — > 77. The uncertainty due to 7r° 
selection is determined from a high purity control sample of J /if} — > 7r + 7r~7r° decays. The 7r° 
selection efficiency is obtained from the change in the 7r° yield in the n + ir~ recoiling mass 
spectrum with or without the 7r° selection requirement. The difference of ir° reconstruction 



efficiency between data and MC simulation gives an uncertainty of 1% per 7r° 20] . The 
uncertainties associated with the EMC energy and \xdE/dx\ requirements are estimated to 
be 0.2% using the control sample J /if) — > 7r + 7r~7r°. According to the MC simulations, the 
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trigger efficiencies for the decays 77/77' '/r] c — > tt + tt~ are almost 100%, and the uncertainties 
are neglected. The systematic uncertainties due to the trigger efficiencies in neutral channels 
v/v'/Vc ~~ >* 7r°7r° are estimated to be < 0.1%, based on cross checks using different trigger 
conditions. The background uncertainties are evaluated by changing the background fitting 
function and the fitting r ang e. Errors on the branching fractions of the J/ip — )■ , yr)/rj'/r) c 
are taken from the PDG 14|. The uncertainty due to the resonance parameters, especially 
the 1] C width, is estimated by varying the parameters by la. Finally the uncertainty on 
the number of J/ip events is 1.3% [7j. Assuming that all of these systematic error sources 
are independent, we add them in quadrature to obtain total systematic errors as shown in 
TableE 



TABLE I: Relative systematic errors (%) of the decay branching fractions. 



Source 




rj/rj' /rj c — > 7r°7r° 


Tracking 


4.0 




Kinematic fit 


2.0 




Photon efficiency 


1.0 


5.0 


7T° selection 




2.0 


MC statistics 


1.4 


1.4 


EMC energy and \xdE/dx\ cuts 


0.2 




Trigger efficiency 




0.1 


Background shape 


4.2/6.3/6.6 


5.6/5.5/8.8 


Branching fractions 


3.1/2.9/24 


3.1/2.9/24 


Resonance parameter 


-/-/8.0 


7-/9.5 


Number of J/ip events 


1.3 


1.3 


Sum in quadrature 


7.3/8.6/27 


8.6/8.5/28 
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VI. SUMMARY 



Since no evident 77, 77' or rj c signal is observed in any decay mode, we determine upper 
limits on the branching fractions B(R — > nn) (R = 77, 7/, r/ c ) with the following formula 

N UP 



N JH , x e x B(J/i; -> 7j R) 
where -A^j 3 is the upper limit on the number of observed events for the signal, e is the 
detection efficiency obtained from MC simulation, where the radiative Ml photon is dis- 
tributed according to a 1 + cos 2 6* distribution and multihadronic decays of 77/7/ '/rj c — > iyty 
are simulated using phase space distributions, Nj/^ is the total number of J/ip events, 
(225.2 ± 2.8) x 10 6 [7|, which is obtained from inclusive hadronic decays. For the decay 
R —7- 7r 7T°, this should be divided by the square of the 7r° — > 77 branching fraction, when 
calculating the branching fraction of R — > 7r°7r°. 

Table Hi] lists the final results for the upper limits on the branching fractions of all the 
processes studied, together with the upper limits on the numbers of signal events, their 
detection efficiencies and signal significances. Here in order to calculate conservative upper 
limits on these branching fractions, the efficiencies have been lowered by a factor of 1 — cr sys . 
For comparison, we also list the best upper limits on r]/r]'/r] c decays to 7r + 7r~ and 7r°7r° 
states to date from the PDG Q. Except for B(r] ->■ tttt), where the KLOE and GAMS-4vr 
Collaborations have huge 77 samples providing the lowest upper limits for 7/ — > 7t + tt~ and 
77 — > 7r°7T , our upper limits are the lowest. These results provide experimental limits for 
theoretical models predicting how much CP and P violation there may be in 77' and 77 c meson 
decays. 
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TABLE II: Summary of the limits on rj/rj' /rj c decays to 7r + 7r - and 7r°7r° states. is the upper 

limit on the number of signal events, e is the efficiency, cT sys is the total systematic error, S is the 
number of statistical significance, /3 UP is the upper limit at the 90% C.L. on the decay branching 
fraction of rjjj]' /t] c to tt + tt~ or 7r°7r°, and B^qq is the upper limit on the decay branching fraction 



from PDG 



Process 


iV UP 

Slg 


£(%) 


CTsys (%) 


S 






R UP 
°PDG 




rj — > TT + 1T~ 


48 


54.28 


7.3 


0.8a 


3.9 x 10 


-4 


1.3 x 10 


-5 


7/ — > 7T + 7T~ 


32 


53.81 


8.6 


O.Ict 


5.5 x 10" 


-5 


2.9 x 10 


-3 


T) c — > TT + 1T~ 


92 


25.27 


27 


1.5a 


1.3 x 10 


-4 


6 x 10" 


1 


7] — > 7T 7r 


36 


23.75 


8.6 


0.6cj 


6.9 x 10" 


-4 


3.5 x 10 


-4 


rf -»■ A 


110 


23.18 


8.5 


2.6a 


4.5 x 10- 


-4 


9 x 10- 


4 


ry c — > ir°ir° 


40 


35.70 


28 


O.lo- 


4.2 x 10" 


-5 


4 x 10" 


4 
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